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Abstract 

In this paper we investigate continuity properties of first and second order shape derivatives of func¬ 
tionals depending on second order elliptic PDE’s around nonsmooth domains, essentially either Lipschitz 
or convex, or satisfying a uniform exterior ball condition. We prove rather sharp continuity results for 
these shape derivatives with respect to Sobolev norms of the boundary-traces of the displacements. With 
respect to previous results of this kind, the approach is quite different and is valid in any dimension N > 2. 
It is based on sharp regularity results for Poisson-type equations in such nonsmooth domains. We also 
enlarge the class of functionals and PDEs for which these estimates apply. Applications are given to qual¬ 
itative properties of shape optimization problems under convexity constraints for the variable domains or 
their complement. 
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Keywords: Shape derivative, shape optimization, convexity constraint, energy functional, Sobolev esti¬ 
mates, optimality conditions, regularity. 


1 Introduction 

In this paper, we focus on regularity estimates for first and second order shape derivatives around non¬ 
smooth subsets of (N >2) for energy functionals involving classical elliptic partial differential equa¬ 
tions (PDE). For instance, we address the following question: given a bounded Lipschitz or convex subset 
Uq C M'^, we wonder what is the optimal regularity of the shape derivatives 

where E' {^Iq) , E" (Qq) respectively denote the first and second shape derivatives around Qq of the shape 
functional i—)• E{Q) = K{x, Ufi, VUn)dx, K{x, •, •) a quadratic polynomial and Uq = Uq{x) the solution 
of an elliptic PDE set in D (see Sections 12.1112.21 for the precise definitions). 

This question, interesting for itself, is in particular motivated by the qualitative analysis of shape opti¬ 
mization problems of the form 

min{J(D), Q C convex, Q G Oad}, (1) 

where Oad is a set of admissible subsets of and J : Oad —)• M is a shape functional which itself involves 
shape functionals D i—>• E{Q) of the above type. 
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The following 2-dimensional example was considered in |19j among other cases: 

J(n) = R{E{Q), |n|) — -P(n), Oad = C open and dQ C {x, a < |x| < 6}| , (2) 

where —)• M is a smooth function, E{0.) is a shape functional related to a PDE, |n| is the Lebesgue 

measure of fl, -P(n) its perimeter and (a, 6) G [0, oo]^. It was proved (see [191 Theorem 2.9, Theorem 
2.12 and Corollary 2.13] and also Section [4] in the present paper) that if the second order shape derivative 
^ !->■ E"{Qo){^,^) around any bounded convex subset Qq is continuous with respect to a norm strictly weaker 
than H^{dU,Q), then optimal shapes of ([I]) are polygonal in {x, a < jxj < b}. In [19], the authors prove that 
such a continuity holds in the two specific examples: when the functional E(i}) is the Dirichlet energy of 
n -that is K = K{x,U,q) = ||g|p — 2/(x) and Uci is solution of the associated Dirichlet problem, or when 
E{Q) is the first eigenvalue of the Dirichlet-Laplacian on D. More precisely, it is proved in [T9| that the 
second order shape derivative of E{-) is, in these two examples and around any open convex domain Dq, 
continuous for the nL°°(9Do) topology (and therefore continuous for the fl^/^'’'^(clDo)-topology 

for any e > 0). The proof of this continuity strongly relies on the 2-dimensional environment. As explained 
below, we prove in this paper that even the full {dQo)-eontinuity holds in this case and even in any 
dimension (see iii) in Corollary 12.6p . Note that this continuity is optimal since, for regular convex domains 
Do, if for instance / = 0 on SDq, then i?"(Do) satisfies 

E"{no){f,^)>cml^/2^anor ( 3 ) 

for any displacement ^ which is orthogonal to SDq. This estimate ([3|) may be obtained by using (5.101) and 
Section 5.9.6. in la¬ 
in Problem ([ll), only convex domains are involved. However, it is also interesting to consider shape 
optimization problems where the PDE is set in an ’’exterior domain” like 

min{ J(M'^ \ D), D C convex, id G Oad}, (4) 

where Oad is as before, but now J involves E{il) = K{Uq,'VUq) where I/q is solution of an elliptic 

PDE set on the exterior domain \ D, which is here the complement of a convex set. It is well-known that 
solutions of such PDE’s are not so regular as in convex domains. These exterior domains are nevertheless 
Lipschitz domains and this is one main reason why it is interesting to look at the case of Lipschitz domains 
even for shape optimization problems involving convex domains like (|4|1. 

In this paper, we use a different strategy to estimate shape derivatives, and we generalize and improve 
the above-mentioned shape derivative estimates in the following directions. 

• A most important generalization concerns the dimension of shapes. In m, the result and the strategy 
were restricted to planar shapes. Thanks to our new strategy, we are able to provide estimates of 
first and second order shape derivatives in any dimension. This is an important breakthrough for 
the study of problems like ([T|) in dimension 3 or higher. 

• We generalize the class of PDE energy-functionals in several ways (see Section 12.11 for precise defini¬ 
tions): the underlying elliptic operator is now a general linear elliptic operator with variable coeffi¬ 
cients, the energy is any quadratic integral functional of the state function and its gradient (and in 
particular does not need to be the energy associated to the PDE defining the state function) and more 
importantly, we consider both interior and exterior PDEs. This last point motivates the next item. 

• We generalize the classes of shapes we consider. Indeed we do not only consider the class of convex 
domains, but we investigate two wider classes, namely the class of Lipschitz domains and the class 
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of Lipschitz domains satisfying a uniform exterior ball condition (we refer to these domains as semi- 
convex domains, see Definition l2.ip . Even if we are interested in applications about convex domains, as 
we pointed out in the previous point, we are interested in PDEs defined in the exterior of some convex 
domain. Whence the consideration of the above kinds of boundary regularity. Estimates for semi- 
convex domains are the same as those for convex domains. They are weaker for Lipschitz domains, 
but probably rather sharp and they are strong enough to be used for our applications (see Section H]) 
which are interesting for themselves. 

• Even in the particular case where E is exactly the Dirichlet energy and 12 is a 2-dimensional convex 
domain, we improve the result of [l9] in obtaining that is continuous in the Lf^/^(c2n)-norm, 

instead of As noticed above, this result is sharp since one cannot expect continuity 

in an iL^-norm for s < 1/2 (even if 12 is smooth). This specific result is actually valid as soon as the 
energy is the one associated to the PDE defining the state function, and is valid in any dimension and 
for any semi-convex domain. 

The method for proving the shape derivative estimates is new. In order to study the derivative at a set 
12o, we first prove adequate estimates of the “material derivative” Ug, where 6 : —>■ is a smooth 

vector held, ' denotes the derivative with respect to 9, Ug = o {I 9), and Ug is the state solution of the 
related PDE in the domain 120 = (/-|-0)(12o). For our purpose, it appears that it is much more efficient than 
using the usual shape derivative Ug. We use the well-known property (see for example m) that without any 
regularity of 12, the map 9 Ug is in general differentiable (even C°° if the involved coefficients are smooth) 
when seen as a ILq (12o)-valued function (whereas 9 Ug is differentiable only when seen as an L^-valued 
function). We show how the regularity of the shape derivatives essentially depends on the regularity of the 
state function Uq. Thus we shall use some sharp regularity results for the solution of a linear second order 
PDE in a Lipschitz or semi-convex domain, in particular VE^’^-regularity results when the data are in 
(see Propositions 13.31 and I3.5jl . 

Note that we insisted in this introduction on the second order shape derivative, but we also obtain 
estimates of the first order shape derivative which seem to be new as well, in this non-smooth setting. 

As an application of these estimates, using the strategy of [THl [19], we prove that any solution 12o of ([I|), 
([4|) is polygonal, if = 2 and if E is in one of two classes described above (where E depends on the solution 
of a PDE in the interior or the exterior of a convex domain, see more precisely Section [2.Ijl . Also, in higher 
dimension, we use our estimates to analyze solutions of the Wdimensional version of Problem (l2|). We 
obtain very strong qualitative properties of optimal shapes, namely that the space of deformations which 
leave it convex is actually of finite dimension (see Theorem 14.Sp . As an easy consequence, we obtain that 
any optimal shape has zero Gauss curvature on any open set where the boundary is smooth. Actually, this 
’’finite dimension” property does contain quite more geometrical information. 

This paper is structured as follows. In Section [2] we introduce the notations, the classes of PDF’s under 
consideration and state our main results. In Section [3] we prove the estimates of first and second order shape 
derivatives. In the last section, we apply these estimates for the analysis of optimal planar convex domains, 
which happen to be polygons and we conclude by analyzing the consequence of our estimates for higher 
dimensional optimal convex domains. 

2 Main results 

2.1 Class of energy functionals 

We deal with energy functionals E(I2) of two main forms. 
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Interior PDE: E{i}) depends on a PDE in the interior of 0, namely 


E{n)= / K{x,UnyUn)dx, 


( 5 ) 


where 


N 


N 


K{x,U,q) = M ^ aij{x)qiqj + aooix)U'^ j +'^l3i{x)qi + j{x)U + (5(x), 

\i,j=i J i=i 

with smooth enough coefficients a, /3, 7 , 6, and Uq solution of 

N N 

U G LU :=-Y^ dj{aij{x)diU) + ^ hdiU + cU = f, in 0. 

i,j=l 2=1 


( 6 ) 


( 7 ) 


Regularity of the coefficient will be made precise later, but they will satisfy the following throughout 
the paper: 


N 

3A>0, VC G Vx G A|Cp < ^ a 7 (x)C*Ci < ( 8 ) 

*j=i 

^ihj) G aij = Uji, Oij = Oji, and c > 0. (9) 

Note that the condition Oij = aji is not restrictive as, in the case the matrix ( 07 ) is not symmetric, 
we can consider = ^( 07 - + aji), which is symmetric, and then in ([ 6 ]) take 07 - instead of aij. Note 
also that no ellipticity condition is a priori required on ( 0 : 7 ). Actually, our strategy can handle much 
more general functional K, see Remark 13.151 

• Exterior PDE: For simplicity, we will restrict ourselves to the model example where 

E{n) = [ |v[/o|^ (10) 

where Uq G solves 

-AC/ = /infl":=R^\a (11) 

and is defined (see [25]) as the set of functions of with trace 0 on dQ, where 






U 


(1 + |xp)^/^(ln(2 + |xp))'^2,Ar/2 


+ < 00 > , 

L2(n-) J 


( 12 ) 


where 52 ^n = 1 if iV = 2 and 62 ^n = 0 if A 7 ^ 2 . 
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2.2 Shape derivative estimates 

In the whole paper, we consider 

(£Br = B{0,R), connected, open, bounded, (13) 

where R> 0 (and large). We set 0 := M'^) the Banach space equipped with the usual norm 

V6» E 0, ||6i||i,oo := sup {||6»(x)|| + ||i:>6i(x)||}. 

x&Bn 

Given a shape functional E{-) : Oad 1^ defined on a family Oad of admissible subsets of M'^, we consider 

f:0^M, £{e)-.= E{{i+ e){no)). 

Then, E is said to be shape differentiable of order m E N* at £Iq (resp. of class near £Iq) if and only 
if T is m times Frechet-differentiable at 0 = 0 (resp. if and only if T is m times continuously differentiable 
in a neighborhood of 0 = 0). When it is well-defined, £'{0) and £"{0) are respectively called the first and 
second order shape derivative of E at flo, and can also be denoted by E'{£Iq) and E''{£Iq). Their values 
at ^,7? E 0 are denoted by T'(0)(^) or E'{Q){^), T"(0)(^,ry) or E''{Q,){^,r]). They are linear and bilinear 
continuous forms on 0 respectively and it is well-known that E', E"(Qq) only depend on the trace on 
d£lo of the deformations ^ and rj (see the proof of Theorem 12.3p . 

Definition 2.1 A bounded open subset Q C is said to he semi-convex if it is Lipschitz and satisfies a 
uniform exterior ball condition in the following sense: there exists r > 0 such that for any x E dfl, there 
exists y E with B{y, r) n O = {x}. 

It is known that a domain is semi-convex if it is locally representable as the graph of a semi-convex function 
(see for example |22l Theorem 3.9]), where a function / : G —)• M is said semi-convex on a convex subset 
C of M” if there exists M E M such that x E C e-)- f{x)-\-M\\x\\‘^ is convex. Recall that £1 is said to be 
Lipschitz if it is locally representable as the graph of a Lipschitz function. 

Remark 2.2 It is easy to check that, if 12 is semi-convex then for Ri small enough, there exist {6i{Ri),a{Ri)) E 
(0,oo) X (0,1) with limi{^_>.o 5i(Ri) = 0 such that 11 is a (<5i(i?i), iT(i2i), i?i)-quasi-convex domain in the 
sense of m- It follows that n is a (d, a, i?)-quasi-convex domain in the sense of [13] for all (d, cr, R) E 
[(5i(i?i), oo) X (0, A:(i2)(T(Ri)] x (0, i?i], with k{R) > 0, k{Ri) = 1, limK^o^(-?2) = 0. 

On the other hand, given a Lipschitz matrix {aifi (i.e. aij E VL^’°°(n) Vf,j = 1,...,A^), then for all 
i?2 £ (0,1), there exists S2{R2) £ (0,oo) with limi^2^.o<^2(.R2) = 0 such that T is a (52(.R2),.R 2 )-vanishing 
matrix in the sense of [13]. It follows that A is (<5, R)-vanishing matrix in the sense of [13] for all (6,R) E 
[<52(i?2),00) X (0,^2]. 

Therefore, we will be able to apply the IT^’^-regularity results proved in Theorem 1.1 of m for the 
solutions to ([7|) in quasi-convex domains (see Proposition 13.5p . Indeed, given p E (l,oo), we first choose 
Ri and R2 small enough so that 6i{Ri), 62{R2) < <5(lV,p) as defined in Theorem 1.1 of [13], and is a 
(<5i(i?i), cr(i?i), Ri)-quasi-convex and (a^) is ((I(i? 2 ), .R 2 )-vanishing. It follows that if R = min{Ri,R 2 } and 
a = k{R)a{Ri) then 12 is a ((5(A^,p), cr, i?)-quasi-convex and (a^) is ((5(A^,p), i?)-vanishing. 

Let us now present the main results of this paper. 

Theorem 2.3 (interior) Let Hq be as in m- For 9 £ Q, we denote by Ue the solution of & in Lid = 

(I -|- 0)(no) (see Provosition \3.tAi and £{9) = E{Q0), where E is given by The following holds. 


5 



i) If atj, bi,c, /, aij, I3i,''y, 6 G m G N*, then [9 i—)• £{6)] is of class near 0 = 0 G 0. 

In the rest of this statement, we assume the previous hypotheses are satisfied for m = 2. 

a) Assume is Lipschitz and L is of the form ©• Then there exists ri = ri(Oo) satisfying ri G (1,2) 
if N = 2 and ri G (1,3) if N > 3 such that, for all r G (ri, oo), there exist Ci = Ci{£lo, L, f, K, r) with 


’■(ano)) (14) 

< C2m\w^ -l/(2r).2r('gQ^), (15) 


for all ^ G 0. 


Hi) Assume Hq is semi-convex and L is of the form ([7|). Then [Iff), (75)) hold for all r G (l,oo). 


iv) Assume Hq is semi-convex, L of the form © with L* = L and E is exactly its associated energy 
(meaning that duK{-, U,VU) — V • dqK{-, U,VU) = LU, for every smooth U). Then [T^ , /IT5l) hold 
for all r G [1, oo). 


Proof. This result is a consequence of Theorem 13.131 which provides an estimate in terms of “interior” 
norms of which after replacing p G (2, pi) by r = p/(p — 2) G (ri := Pi/{pi — 2), oo) gives 


\E'{Il,m\ < Ci||e||u.i,.(no), |i7"(b!o)(C,OI < C^2||C|lW(no)- (16) 

If £Iq is semi-convex then we can have any p G (2, oo) which leads to any r G (1, oo), and if moreover L is 
self-adjoint and K is the associated energy, then we can have p = oo,r = 1. 

Now, we use the well-known fact, that |) depend only on the values of ^ on d£lo. More 

precisely, 

=> E'illom = E\nom, E"ino){^,0=E"miii). 

Indeed, let Cn be a sequence in converging to ^ | in lT^’^''(no)- Let Zn G C'^([0, oo) : be 

the solution of 

Vt > 0, Vx G M^, dtZn{t,x) = (n{Zn{t,x)), Zn{0,x) = X. 

Note that Zn{t, •) G 0 (see [E]), and for x in some neighborhood of Sflo, Zn(t, x) = x for all t > 0. It follows 
that Zn{t,Ilo) = rio, where Zn(t,Q,Q) = {zn(t,x), x G Hq}, and therefore E{zn(t,QQ)) = E{Qq) for all t > 0. 
In particular 

E'ino){^-i) = hm^ E{znit,no)) = 0, 

n^oocit |t^o 

E"m{C-iC-0 = E{zn{t,no)) = 0. 

n^oocLE |t=o 

From this property and (flbll . we deduce 

|i7'(L!o)(OI <Llimf{||?1lH^iv(no) i = iondllo}, 

and the similar property for E"{Qq){^,^) with the lT^’^”-norm. We then apply Lemma l3.II to deduce (I14I) . 

P^. □ 
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Remark 2.4 The result i) of Theorem 12.31 a priori implies that E'{Qq) and are continuous on 0 

for its iT^’°°(no)-norm (as a linear and a bilinear form respectively). Using the previous trace analysis, 
we deduce that they are also continuous for the lU^’°°((90o)-norm. The estimates (I14p . (jl5|] improve this 
a priori information and show that ^^'(^ 0 )) are continnons in and 

respectively: note that the regularity gets stronger as r gets smaller. 

In case iv), the estimate (I15p is valid with r = 1 (which yields the R^/^-continuity) and cannot be improved 
since, as already noticed in the introduction (see ([3])), E"(i}Q) has i4^/^-coercivity properties. 

Again in case iv), it is interesting to notice that (fHI) with r = 1 can be stated as U'(flo) € (L^((9Qo))^ = 
L°°{dQo). Again this is sharp in general; if for example E is the Dirichlet energy, 

E{n) = l [ \VUn\‘^- [ fUn, 

^ Jn Jn 

it is well-known that E'{^q) = which indeed belongs to L°°{d^o). 

Theorem 2.5 (exterior) Let Hq be as in and Lipschitz and let f G IU™’°°(M^), m G N*, with 
supp(/) d Bji. For 9 £ Q, let Ug G HQ^{Llg) be the solution of /flil) in (see Provosition \3.16\) and 
let £{6) = E{Llg) with E given by HU) . Then £ is of class near 0 = 0 G 0. Moreover, there exists 
ri = ri(flo) satisfying ri G (1,2) if N = 2 and Ti G (1,3) if N > 3 such that, for any r G (ri,oo), there 
exist Ci = Ci{£lQ, Uq, R, /, r) with 

\E'{nom\ < CiiiciIvKi-id, ^{dno) */ 

\E"{noM,0\ < C2mli -l/(2r),2r(9Qp) if 171 >2, 

for any G 0. 

Proof. Similarly to the proof of Theorem 12.31 this result is a consequence of Theorem 13. 191 combined with 
the trace result of Lemma EH □ 

Previous estimates can be written using a different range of Sobolev spaces. In the following statement, 
we write these estimates in the R^-norms, which are relevant for our applications (see Section E]). 

Corollary 2.6 Let Uq be as in m- 

i) Under the same conditions as in ii) of Theorem \2.tA or as in Theorem \2 . 51 and with ri = ri(no) os 
in these theorems, we have 

|.U"(Llo)(C,OI < <^2 ||^|||/s(90o), Vs G (S2,l], = min|l, , ( 20 ) 

for every ^ G 0. 

ii) Under the same conditions as in Hi) of Theorem \2.3l ifTPI) (resp. holds for every s G (0,1] (resp. 

s£ [1/2,1]). 

Hi) Under the same conditions as in iv) of Theorem \2.3l UM (resp. flOj) ) holds also for s = 0 (resp. 
s = 1/2). 


I-E'^(^o)(0l < Vs G (si, 1], 


Si = mm 


(17) 

(18) 
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Proof. Inequalities (fT9l) . (1211 follow from (fT^ and (fT7|l respectively and from the embedding result (f2T]l 
recalled below which may be obtained from [301 Th. 1.107]. Indeed, OQq is an {N — l)-dimensional Lipschitz 
manifold and since we deal with exponents less than 1, these embeddings carry over to Lipschitz manifolds. 


Vti,t 2 G [0,1], V(pi,p 2 ) € [l,oo], 

r N -1 N -1 

ti - > t 2 - and ti > t 2 

I Pi P2 

We apply this result with ti = s E [0, l],pi = 2 and, [t 2 = 1 — l/r,p 2 = r] to obtain ([T9]l and [t 2 = 
1 — l/( 2 r),p 2 = 2r] to obtain ([20]) . 

For the proof of ii) and iii) we apply the same embeddings and iii), iv) from Theorem 12.31 □ 

Remark 2.7 The restrictions on si, S 2 in the Lipschitz case i) of Corollary 12.61 are coming from the restric¬ 
tion p < Pi = Pi (Ho) as indicated in Theorems 13.9113.181 We do not know the exact value of this Pi, but 
we know that pi > 4 if = 2 and pi > 3 if > 3. Thus (recall that ri = Pi/(pi — 2)), we may write 
ri=2 — fTifAi = 2 and ri = 3 — fi if A^ > 3, for some small cr > 0. 

In the range of L7®-spaces, the best estimates we obtain in the case i) of Corollary 12.61 mav also be written 
as follows where the constants Ci,C 2 here do not depend on ^ E 0, but may depend on other variables, 
especially e): 

• if A^ = 2, then there exists e > 0 small such that 

(22) 

|i?"(^2o)(C,OI <C2||C|lii-.(9no)- (23) 

• If A^ = 3, then there exists e > 0 small such that 

|i^'(^2o)(0l <C'i||^||^^i-.(aoo)- 

3 Estimates of shape derivatives 

3.1 Description of the method 

We describe in this paragraph the main idea of the shape derivative estimates on the model example 

E{n)= [ K{Un,VUn) := [ hvUn\\ 

Jn Jn ^ 

where Uq is solution of d?]) with L = — A and II, Oq are as in (|13l) . The main point of our approach is that we 
will make interior estimates involving VF^’'^(IIo)‘iiorms of the directions ^ E 0 of differentiation. As already 
explained in the proof of Theorem 12.31 we will then use the trace Lemma 13.11 to obtain estimates in terms 
of VF^“^/'^’'^(dno)-norms of By doing so, we avoid integrations by parts which would be impossible due 
to the poor regularity of the boundary of Hq. Moreover, as we will se below, the shape derivative estimates 
will only depend on the regularity of I/q := 

One first has to show that E is differentiable. 

By changing variable x = {I + 0){y), for 0 E 0 small, we have 

^(^)= f l^Ue-Me-VUedx, 

Jno 2 


W^^’P^idno) C W^^’P^dflo). (21) 








where 


Ue = Ueo{I + e), Me = [l + ve]-^\l + ve]-^det[l + vei Ue = Une- (24) 

Note that 9 i—)• Mg is C°° from 0 to near 0 = 0, see [HI |28]. Then the differentiability of 

9 1 -^ £(9) fully depends on the differentiability of Ug. As it is classical, under reasonable assumptions on L 
and the data, by using the implicit function theorem, one can prove that 9 £ Q Ug G Hq{Qo) is of class 
C™ near 0 = 0, which implies that £ is of class C™" near 0 = 0. 


Next, let us fix a direction ^ G 0; we will denote by (•)' the derivatives with respect to 0 in the direction 
In general, Uq is more regular than H^{Qq) and we consider p G (2,oo) such that Uq G (see 

Propositions 13.31 and 13.51 below). We may write: 

£'(0m= f VUg-Mg-VU'g + l-VUg-M'g-VUg, (25) 

jQo ^ 

which at 0 = 0 implies fsee ffMlP 


i^'(o)(e)i < 

< 


c 

c 



Vf/o • VC/^ 
Vf/o • VC/' 


+ [ |VC/o|2|VC|) 

+ ll^o||^i,p(OQ)||^||pvup/(p-2)(no) 


(26) 


% 


The second term in the last estimate being satisfactory, it remains to estimate the first term which involves 
Note first that it would not help much to use a too simple Holder inequality like 



VC/o • VC/' 


< C'II^o||h1(Oo)II?IIiV1'°°(Oo)’ 


which uses only the starting information that C/q is a linear continuous map from 0 to 

Nor is it appropriate to write the term C/q in terms of the shape derivative C/q, i.e. C/q = C/q + VC/q • 
Indeed, in such a case one would have 



VC/o • ^Uo 



VC/o • VC/^ + VC/o • D^Uo • -e + VC/o • V^ VC/o. 


But we would need here regularity for D^Uq and it is not available for the case we are interested in. Even 
if in the semi-convex situation, we can get some significant information on the first derivative, it becomes 
quite more difficult for the second derivative (see however m for some progress in this direction). 

For these reasons, we proceed with the estimate of the term with C/q by going back the state equation : 
LUe = /) Ug £ HQ^Qg). Recall that we chose L = —A for simplicity here (but the ideas will be the same 
for general L). The weak form of LUg = f in £lg transported in Hq is 


/ 

J Qq 


ipLgUgdx = / fgipdx, for all ip £ Hq{£Io), 


iQq 


with fg = fo(I + 9) det[/ -|- V0] and 

Lg : ^ H-\no), Lgcj) = V- {Mg • VcP), 


(27) 


(28) 


where Mg is defined in (12^ . Note that Lq = L and {LgUg)' = LgJJ'g + LgUg. By differentiating (1771) with 
respect to 0 in the direction we obtain 




pLgU'dx = / pif's-L'gtjg), 

J r^o 


(29) 
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(30) 


which at 0 = 0 gives 

[ = f ^(/'-V-(M^V[/o)), 

J Qq J Qq 

The estimate (fM]l suggests to take (p = Uq £ in pOl) . which then yields 



V[/o • Vi7' 


< Cinoj) [ (|c/ol + |vc/op)(|c| + |vei) 

J Qo 

< C{Q.q, f ,p){l + ||C^o||vKi-P(no))ll^o||vKi’P(r2o)ll?lllVi’P/(p-2)(Oo)- 


Thus, going back to we deduce 


(31) 


|T'(0)(OI < ^^(f^o,C/o,/,p)||ellTvW(P-2)(no)- (32) 

For more general K and L, the choice of p is more involved (see the proof of Theorem l3.9p . but the procedure 
is the same. 

In the same spirit, that is to say by differentiating (12511 and (12911 at 0 = 0, we can obtain a similar estimate 
for |T"(0)(^,^)|. For more details see Sections lOl 13.41 


3.2 The trace lemma 

Lemma 3.1 Let he as in l[TM) and Lipschitz. Let also q E [l,oo],s E 


T 1 

,1H— 1, or q = s = 1. Then 

there exists C = C{s,q,Ll) such that, for every ^ E IT^’°°(fI, M^) n IT^’'?(n, M^), 


mf{ miws.in), i = ^ on dLl] < (33) 

The main tool for the proof of this lemma is the following classical trace/extension Theorem: 


Theorem 3.2 Let LI he a set with Lipschitz boundary, q E [1, oo] and •s E (|, 1 + |) or q = s = 1. The trace 


mapping Tr initially defined on C°°{Ll) extends as a hounded linear operator from W‘^’^{Ll) to 
Moreover there exists a bounded linear operator 

Ext : W"~T'^{dn) 

such that Tr o Ext = Id on W^’'^{dLl). 

See |12l Theorem 3.1] and the references therein for the case (7 = s = 1, and for example |14l Thl pl97] for 
the other cases. 

Proof of Lemma 13.11 Let ^ E IT^’°°(n, M'^)nIT^’'^(n, M^) (we drop the notation in the following). We 
remark that E VF* Using the extension defined in Theorem l3.2[ we can define ^ := Ext(Tr ^) E 

VF^’'^(fl) satisfying 




(34) 


where C = C{s,q,Il). 

Therefore ^ ^ E IUg’‘^(0) and is therefore the limit in of functions E (7^(0). In other words, 

I is the limit in IT®’'^(n) of pn = f + cun- Clearly pn £ IU^’°°(II) n W^’'^{Ll) and Tv{ipn) = Tr(0) so that 

inf I |||||u/'>.9(n), I E IU^’°°(II) n | ^ on 5II | < ||(/::n||rv'»,9(o), 

and letting n ^ 00 and using (j34ll . we obtain the first estimate (l33ll . □ 
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3.3 Shape derivative estimates for an interior PDE 

In this section we will prove estimates for the first and second order shape derivatives of the energy (j5]) . 
As explained above, it will rely on estimates of Uq. 

3.3.1 Existence, uniqueness and regularity of the solution Uq 

Proposition 3.3 Let be as in m and Lipsehitz. Let L be as in & with the matrix (aij) satisfying 

mM- 

i) Assume bi,c E L°°{Ll) and c > 0. /// E then there exists a unique solution Uq E Hq{LI) of 

&■ 

a) Assume aij E bi,c E L°°{Q) and c > 0. Then there exists pi = pi{Ll) satisfying [pi > 4 

if N = 2 and Pi > 3 if N > 3], sueh that for every p E {p'i,pi) (where 1/pi + = Ij o,nd every 

f E W~^'P{Ll) the problem 0 admits a unique solution Uq E IPQ^’^(n) and 

||%||ivi.p(n) < C{Ll)\\f\\iY-i,p(^Q). (35) 

Proof. The point i) is standard, see [8]. The sharp regularity result in ii) has been first proved for 

L = —A in m Thm. 0.5]. The complete proof of ii) is based on [29] Theorem C] and Bemark 13.61 for 

Lu = Remark 13.71 for L of general form. See also Remark 13.81 for a different proof 

with a slightly stronger regularity on Uij. Note that, according to [29], the same result is valid with quite 
weaker regularity on Uij, like asking that they be in VMO(M'^). □ 

Remark 3.4 Notation: If L satisfies the assumptions of ii) above, then so does its adjoint L* as we can 
easily see by writing {biu)xi = {bi)xiU + huxi, where bi E IT^’°°(n). Let pi, p'^ (resp. pi, p'f) be the 
numbers associated to the equation LU = f, U & Wq’^{Q) (resp. L*U = f, U ^ Wg^’^(II)) as given by ii). 
Proposition 13.31 Then we set 

Pi := mm{pi,pl}, p[ := max{p^,p']]}. 

Note that, if Ll is Lipsehitz, then p;^ > 4 if A = 2 and p^ > 3 if A > 3. 

Proposition 3.5 Let II be as in m and semi-eonvex. Let L be as in & with the matrix (aij) satisfying 

mM- 

i) Assume aij E IT^’°°(n), bi, c E L°°{Q) and c > 0. Then for every p E (1, oo) and for all f E 
the problem admits a unique solution Uq E ITQ^’^(n) and (E^) holds. 

ii) Assume aij E IL^’°°(0), bi,c E L°°{Q) and c > 0. Then for every f E L°°(n) there exists a unique 
solution Uq E Hq{Q) n IT^’°°(n) of ([^. 

Proof. The result in i) was proven for L = —A, see (6] Corollary 1 and Remark]. For general L, we 
refer to m as explained in Remark 12.21 and we use Remark 13.71 below to include the first order terms. For 
ii), as / E L°°{Ll), from i) it follows that Uq E VF^’^(n) for all p E (l,oo), which implies Uq E 
Furthermore, applying [SI Thm. 8.8], it follows u E Applying [S] Theorem 9.13] (locally in Q, with 

T = tj)) and then (8] Lemma 9.16], we get u E for all q > 2. From Sobolev embeddings it follows 

u E (^^(n) n Then ii) follows from [8l Thm. 15.9]. □ 


11 


Remark 3.6 Theorem C in [29j states that if Uij G VMO(M'^) (and in particular if aij G 
and L = —di{aij{x)dj), then T := is continuous from LP{Q) to {LP{Q))^, for p G (l,Pi), with 

Pi = Pi(2) > 4, pi{N) > 3 for N > 3. Then it follows that L is continuous invertible from ITQ’^(n) 

onto for p G {p'i,pi)- Indeed, first note that this is equivalent to show that V(L)“^div is continuous 

from to Next, we note that from [29l Theorem C], T* = (L)“^/^div is continuous from 

(L'^(n))^ to L'^(n), for q G (y]^,oo). Then the desired continuity for V(L)“^div follows from the fact that 
V(L)-idiv = TTT 

Remark 3.7 To complete the proof of ii) in Proposition 13.31 and of i) in Proposition 13.51 . let us check that 
if the map u —Au := —di{aij{x)dju) defines an isomorphism from ITQ^’^(n) into then so does L. 

Indeed, \ei B = L — A, Bu = hdiU + cu, and consider the equation 

Lu = f, uGW^’^in), few-^’P{n), (36) 

which is equivalent to Au + Bu = f. Multiplying this equation by A~^, we get the equivalent equation 

iI + K)u = g, g = A-^feW^'Pin), ueW^'^in), (37) 

where K = A~^B : ITQ’^(n) i—ITg’^(fl) is compact since B is compact from ITg’^ into W~^’^ and A~^ is 
continuous from into IPo’^- 

Then I -\- K satisfies the Fredholm alternative (see e.g. [H Thm. 6 . 6 ]). Furthermore, Ker{I + K) = 
Ker{L) = {0}. Indeed, assume u G Wq'^{Vi), Lu = 0. If p > 2, we can directly use the uniqueness result 
in [U Thm. 8.1 and Cor. 8.2] to deduce u = 0. If p G (1,2), we obtain the same conclusion by using that 
L* : I—>■ is onto which follows from [ 8 l Thm. 8 .3 and 8 . 6 ] (see also the remark following Corollary 

8.7 in [ 8 ]). 

Then, Fredholm alternative [H Thm. 6 . 6 ] implies that I K defines an isomorphism from onto 

itself. □ 

Remark 3.8 Different approaches may be quoted for the last point of Proposition 13.31 They only require 
a slightly stronger regularity hypothesis (like aij G C'^''"^(D), 7 G (0,1)), but they are also quite interesting. 
Let us just recall the idea: we need to focus only in the case L = —di{aij{x)di), because the case of gen¬ 
eral L follows from Remark 13.71 The case N = 2 follows from [2T] . For the case N > 3 we consider the 
manifold {Q,g), where g is the metric given by <7 = det(aij)^/*'^“^^^“^. Then (from |24( p. 186] and the 
formula of Laplace-Beltrami operator, see [2ll (1.1)] for instance), the equation Lu = f, u ^ Wq’^{Vi) with 
/ G VF“^’P(D), becomes —A^gU = h, with h = det(ajj)^/(^“'^)/ G VF“^’P(D). Note that g G ^nd 

VF“^’^(D) norms of / and h are equivalent. Then the claim ii) follows from [23( Corollary 13.2]. 

3.3.2 Estimates of the shape derivatives of the solution 

The results of this section give estimates for Uq and Uq. Note that while for Uq we prove an 77^ (Dq)- 
estimate, we also prove more involved estimates for Uq and Uq. These estimates are motivated by the 
functional K involved in the energy ([5]), see also Theorem 13.131 

Theorem 3.9 Let Qq be as in I113\} and let L be as in & with the matrix (uy) satisfying Let 6 £ Q, 

Qq = {I + 0)(Do), Uq G Hq{Q.q) solution of & (see i) in Provosition \3.3\) . Uq = LfQO {I + 6). 

i) If aij,bi,c, f G W'^'°°{Br), m G N*, then JJq G C'”*(V, 77q (Dq)) where V is a neighborhood of 9 = 0 in 

0 . 
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In the following, we assume the previous hypotheses are satisfied for m = 2. 


a) If moreover r^o is Lipschitz and aij E /3j,7,Q;oo 


/ duK{;Uo,VUo)U[, + dgK{;Uo,VUo)-VU[, < 


/ Qq 


duK{; Uo,VUo)U(f + dgK{;Uo,VUo) • vu;f 


< 


E L°°{Q,) then for all p E {2,Pi) we have 


Cl ? VF1’2p/(p-2)(Oq), 

(38) 

Cl\\i\\w^’P/(p-^) (Ho)’ 

(39) 

C 2 11 11 1,2p / (p - 2) ( J-J JJ p 

(40) 


where Pi is defined in Remark \3^ (pi > A if N = 2, p^ > 3 if N > 3), Uq = Ug{0){f), Uq = (0)(^, ^), 

Ci = Ci{no,LJ,K,p). 


Proof. 

Step 1. The proof of i) is classical. It is based on the implicit function theorem. 

After changing the variable x = (/ + 0){y), the weak form of ([7|) in klg is transformed in Qq as follows 


{Lg(le,ip)H-i{no)xH^{no) — / E flo(no) 


Ho 


(41) 


where Lq : HQ{klQ) i—)■ H ^(flo) is defined by 

{LeU,ip) = {di{mdkU) + bkdkU + cU,^) 

akidkUdiif + {bkdkU + cU)ip, for {U,ip) E Fo(^o)^, 


' Ho 


(42) 


and 


^ki = 'EijO'ij ° + ^)TkiTijdet[I+ V9], bk = 6^ o (/+ 6»)rfcidet[/+ V6»], 

c = CO (/ + 0)det[/ +V0], / = /o (/ + 0)det[I +V0], 

T = {Tki) = [id + ve]-\ 


Consider F defined by 

F: QxHli^o) 

{e,u) 


(f {LelJ - f,ip} 


It is easy to check that F is well defined and of class C"* in a neighborhood of (0, Uq) E 0 x Fg (Uq), and 


duF{0, Uq){U) : ip {LU, p)= [ + (bAU + cU)p, Vt/, ip E Fg^Ho). 

J Qq 

Note that from Proposition 13.31 rIf/F(0,[/o) defines an isomorphism from Hq{'D.q) to H~^{'D.q). Then, from 
implicit function Theorem there exists a map, U : 6 U{9), such that F[9,U{9)) = 0 for ||0||0 small. 

From the uniqueness of solution to ([7]) we obtain Ug = IJ{9), which proves the regularity oi 9 ^tJe. 

Step 2. We prove the estimates (f38]l - (f40]l by differentiating (j4T]l with respect to 9. The differentiation is 
allowed because Ug is differentiable and aij,bi,c,f are regular enough. 

Differentiating (I41h once with respect to 9 gives 

{LglJ'^p) = {f'-L',Ug,p). (43) 
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We take 0 = 0 in (j43|) . Note that Lq = L. By isolating all the terms with Uq we obtain 

= {^'ki^kUodup + {h'kdkUo + CQUo)ipj . (44) 

We now choose suitable test functions (p to prove (l38]l and (f39]l . Let p G Hq{Qq) be the solution of 

L*p = jj'^-^U'^=:g\ (45) 

Note that g* G H~^{U,q) and then the solution ip is uniquely defined in Hq{VIq) and satisfies (see [T7] l 

\H^{no) ^ < C'IIL'oIIhLW)’ ^ = C{^lo,L). (46) 


Then from (fUll it follows 

\\u(,\\mm < Cl [ (i + |c/o| + |vc/o|)(iei + |v^|)(|(^| + |V(^|) 

J Qq 

< (^ 1(1 + ||[/o||wi-P(r2o))ll'^llvUi'2p/(p-2)(no)||</5||Hi(f^o)’ 

where we applied Holder inequality with three terms, and Ci = Ci{Q-q,L, f), p G (2,pi), pi = pi(Ho) as 
given by Proposition 13.31 for L. Together with (HOl) this proves ([38]) . 

Now we consider p G iLg(Ho) the solution of 

L*p = duK{;Uo,VUo)-V-dgK{;Uo,VUo) 

= iaooUo + 7)-V ■iaVUo + P)=:g*. (48) 

Let q G {2, pi). Recall that Uq G W^’'^(Ho) by Proposition 13.31 From ([6]), we deduce 


duK{; Uo,VUo),dqK{; Uo,VUo) G L''(L!o). 

It follows that g* G IT“^’'?(Ho). By Proposition 13.31 applied to L*, the solution p of (|18]) is uniquely defined 
in ITg^’'^(Ho) and satisfies 


ll<^llwi.'j(Qo) ^ C\\g*\\w-i,q(^Q^) — C, C — C{Q.o,L,K,Uo,q). 


(49) 


Then from (fill) it follows 


[ duK{;Uo,VUo)U;i + dqK{Uo,VUo)-VU'^ 

J 


{L*p,U'^) ={LU',p) 


< Cl [ (l + |17o| + |VHo|)(|^| + |Ve|)(|(/.| + |V¥.|) 

J Qq 

< < 71(1 + ||Ho||vi/i.p(n(,)||(/7||u/i.'j(fio)ll'^llwi-i/(i-i/p-i/9)(r2o))) 


(50) 


with Cl = (71(00, L, K, /, q). Here we chose q '■= p where p is given in ii) of Theorem 13.91 and the estimate 
follows. (Note that the use in considering q ^ p in the previous computations will appear later in the 
case of semi-convex domains). 

Step 3. Differentiating (j43j) at 0 = 0 and isolating the terms with Uq gives 


{LlJ'^.p) 



i^ki^kUodip + {b'k^kUo + CqU(i)p) 



(^b'kidkUodip + {b'i.dkUQ + CqUq)p'^ . 


(51) 
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It implies 


(LV,I7^')| < C2^ ((l + |C/o| + |V[/o|)(|C|2 + |VC|2) + (|C/'| + |V[/^|)(|C| + |Ve|))(|¥.| + |V</p|), (52) 


where C 2 = C 2 (fio, L, /). 

Then ([^ with (p solution of (USD, together with (fHSp yield 


' Ho 


duK{; Uo,VUo)U(^ + dgKi; Uo, VC/o) • VI/^' 


< C'2f (1 + ||C/o||vKi.p(no))lll?l + llLi/(i-i/p-i/9)(no) + ll^^ollL2(f2p)|||||v(/l,l/(l/2-l/9)(OQ) 




IVKi’9(Oo) 

< C2il + \\Uo\\ VFi’P(no))ll^ollvKi.5(no) 

|^l,2/(l-l/p-l/q)(QQ) + ||Cllvyl.2p/(P-2)(no)ll?llvFl’29/(9-2)(no) 
with C 2 = € 2 (^ 0 , L, /, K,p, q). Again we choose q := p and this proves (HUp . 


(53) 

□ 


Proposition 3.10 Besides the assumptions of ii) in Theorem AS.fA we assume Hq is semi-convex. Then the 
following holds. 

i) For all p G (2, 00 ) we have 


\\U()\\h^{Qo) 

[ duK{; [/o, VC/o) 17 ' + d,K{; Uo,VUo) • vO' 

J Qq 

[ duK{; Uo,yUo)U'f + dqK{; Uq, VC/q) • VU[f 
J r^o 


< 

C*! C Lf/pOo)) 

(54) 

< 

^ n/i-p/(p-2)(Oo)’ 

(55) 

< 

^2 ||^||^l,2p/(p-2) (Oq)) 

(56) 


where Ci = Ci{FlQ,L,f, K,p). 

ii) If furthermore L is self-adjoint and E is its energy associated, i.e. duK{-, U, VU) — 'V-dgK{-, U, VC/) = 
LU, then we can take p = 00 in (ED, (ED- 

Proof. We proceed as in Theorem 13.91 using the extra property that flo is semi-convex. We now have 
Uq G W^’°°(r2o) (see Proposition 13.51) . Moreover, the solution tp of (HBIl satisfies (p G for all 

q G ( 2 , 00 ). Then we can take p = 00 in (f47ll which implies (f54ll . Furthermore, we can apply (|50D and (f53D 
with p = 00 and q G (2, 00 ) arbitrary, which prove ([5^ and ([SOp . 

In the case 12 semi-convex, L self-adjoint and K the energy associated to L, then L = L* and L*ip = LUq. 
Therefore, p = C/q in (Hsp . Hence Uq G H/^’°°(I2o) and p G VF^’°°(no)- Then we proceed as above with 
p = q = 00 . □ 


Remark 3.11 The results of Proposition 13.101 hold for any Uq such that Uq G H/^’°°(I2o) and such that 
p G IFQ^’'^(no) L*Lp G H/“^’'?(I2o) is bounded and invertible for all q G (2, 00 ). 

Remark 3.12 In Theorem 13.91 we could try to estimate ||C/g Indeed, in step 3 of Theorem 13.91 we 

take ip to be the solution of L*ip = Uq — AC/q =: /*, so that ||f^o IIhi(Oo) ~ Uq). Then we can proceed 
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as in (l53]l . However, as we have only Uq,Uo G then f* £ H ^(Oq) and so (/? G Hq{Qq) (sharp in 

general). Therefore, in ([531) we have q = 2 and it leads to 


l^oll^i(f^o) - ^2(1 + ||h^o||vi/i.p(no))ll^o||ivi. 2 (no) 


(57) 


which is not appropriate as it contains the strong norm ||Cllivi’'=°(no)- 


3.3.3 Estimates of the shape derivatives of the energy 

Theorem 3.13 Let flo be as in l[TM) and let L be as in ([^ with the matrix {aij) satisfying The 

following properties hold. 

i) If aij,bi,c, faij,j3i,'y,5 G m>l, then £ : 0 £{9) is of class C™ in a neighborhood of 

0 = 0 in Q. 

In the following, we assume the previous hypotheses are satisfied for m = 2. 
a) If Hq is Lipschitz, then for all p G (2,Pi) we have 


l^^(0)(OI ^ C'l||^||^yl,p/(p- 2 )(f^g), (58) 

\£"{0M,O\ < C 2 ||e||^i, 2 p/(p- 2 )(f,„), (59) 

where Ci = L, f, K,p) and pi (introduced in Remark \3.4^ satisfies pi > 4 if N = 2 and Pi > 3 

ifN>3. 

Hi) If Hq is semi-convex, then £5^) and /f5^) hold for all p G (2, 00 ). 

iv) If rio is semi-convex, L is self-adjoint and K is its associated energy, then and l5^) hold for 
p = 00 . 


Remark 3.14 Theorem 13.131 can be stated in a more general form, based on the regularity of the state 
solution related to the operator L and its adjoint L*. Namely, if 

a) aij,bi,c,f, aij,fii,j,5£W‘^’°°{no), 

b) the operators L and L* define isomorphisms from ITg’^(ll) 1 —>■ for some p G [l,oo] 

then (f58]l and (fM)) hold for this p. 

Proof, [of Theorem 13.13] 

Step 1. Note that by changing the variable y = (/ + 9){x) we have 


£(6) = / K{I+ e,Ue,^[T]-VUe)det[I+ Ve]dx 


f 

Ino 


-^dkidkllediUe + -dooUg + fikdkUe + ^IJe + < 5 ), 


where 


dki = oiij o {I + e)TkiTijdet[I+ V9], fii = fiio {I + 9)Tkidet[I+ V9], 

7 = 7 o (/+ 0)det[/+ V0], (5 = 5 o (/+ 0)det[/+ V0]. 


( 60 ) 
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The differentiability of £ follows from the regularity of K and i), Theorem 13.91 
Step 2. Assume ffo is Lipschitz. Differentiating (I60p gives 

£:'(0)(O = [ akidkU'sdiUe + aooUeUe + MkU's + We 

J Qq 

+ [ (—^ki^kUedilJg-\-—aQQUQP'j^dkUe-\-j'Ue + 6'. (61) 

Jno ^ 

Taking 0 = 0 in (f6T]l and using (1391) gives 


l^'(o)(OI < 


+ 


[ duK{x,Uo,VUo)U() + dgK{x,Uo,VUo) 
J Qq 

Cl [ (i + |C/o|2 + |V[/o|2)(|e| + |VC|) 

J Qq 


VU', 


^ C'i(||^||yi..i,9/(9-2)(Qg) + (1 + l|t^o||^i,p(Oo)ll^llvyi’P/(p-2)(no))> 


(62) 


with Cl = Ci{no,L, f,K,p) and p,q£ (2,Pi). 

Differentiating (16111 at 0 = 0, isolating the terms with Uq, using the lT^’°°-regularity of the coefficients of 
L and K, the fT^’P(Do) regularity of Uq and the estimates (f38]l . (flO]) yields 


l^"(o)(?, 6 l < 


f 2 o 

+C 2 


duK{x, Co, VCo)C^' + dqK{x, Co, VCo) • VC^' 


'^0 - 
|2 


It^or + IVC'I^ + (1 + |Co|^ + |VCo|^)(|e|^ + |vcn 

2 Mltl|2 


^ C2(||?|liyl,2g/(Q'-2)(J-J|j) + (1 + l,2p/(p-2) (llg) 


(63) 


with C 2 = C 2 {£Io,L, f,K,p) and p,q£ (2,pi). 

Taking q = p in ([62]) , (l63|) proves ii). 

Step 3. If Do is semi-convex then from Co G IT^’°°(Do) and i). Proposition I3.10( we can take p = 00 and 
q € ( 2 , 00 ) in (f62]l and ([63]) . which proves hi). 

Step 4- Finally, if Dq is semi-convex, L = L* and K is the energy associated to L then from Co G VF^’°°(Do) 
and ii). Proposition 13.101 we can take g = p = 00 in ([62]) and (l63]l . which proves iv). □ 


Remark 3.15 The technique we have used to obtain estimates for C'(D) and C"(D), can be applied to 
shape functionals E defined as in ([5]) but involving more general K. Once the differentiability of E is proven, 
our computations can be used similarly to get estimates whose exponents will depend on the growth of K 
at infinity. Note that the proof of the differentiability of E may rely in these cases on the differentiability 
of the map 9 £ Q U{6) G VF^’^(Do) with certain p>2, which is not known as far as we know. 

3.4 Shape derivative estimates for an exterior PDE 

In this section we will apply the technique described in Section 13.11 to estimate the first and second order 
shape derivatives of the energy (11011 . related to the problem ([1111 in the exterior of a domain D. 

Even in the case when D is convex, its exterior D® is just a Lipschitz domain. Therefore, subject to the 
existence, uniqueness and regularity of the solution of the problem (jlip . which represents some particularities 
as it is in the exterior of D, the method we developed in Section 13.11 and used for the interior problem in 
Section 13.31 also applies for the problem in the exterior. 

As in the previous section, we consider 0 = and all the domains D under consideration 

will be assumed to satisfy at least m- 
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3.4.1 Existence, uniqueness and regularity 


To analyze problem (fTTI) . it is appropriate to consider the space as introduced in (fl^ . This space 

is a Banach space, see [25], and even an Hilbert space if equipped with the inner product 


{U,V)Huo^n^y.= [ ( 


UV 


Let denote the closure of in It can be shown that equipped with the 


(1 + |xP)(ln(2 + |xp)'^ 2 ,Ar 


rl,0 


+ VU ■W] dx. 


(64) 


1 , 0 / 


inner product 


iU,V) 




-I 


= / VU-VVdx, 


(65) 


is an Hilbert space, see |25] . Furthermore, the norms in generated by the inner products (|64p and 

(1651) are equivalent. Let denote the dual space of 


For the solution of (jlip we have the following regularity result. 

Proposition 3.16 Let Ll be as in m- Assume it is also Lipschitz and simply connected. Let f € L°°{M.^), 
supp(/) C Br. 

i) Then there exists a unique weak solution U E LfQ’°(H®) to /fIT]) satisfying 

||C/kl, 0 (Oe) <C(H,i?)||/||^-l, 0 (ne). ( 66 ) 

a) There exists pi = pi(Ho) satisfying [pi > 4 if N = 2 and pi > 3 if N > 3], such that U E Wq'^{LI^) 
for any p E {p'i,pi) where l/pi + l/p^ = 1, and 

\\U\\wiyny<C{n,R,p)\\f\\L^, n% = n^nBR. (67) 

Proof. 

Step 1. Note that the weak solution U of (fTT]l satisfies 

[ VU-V^= [ fif, Vy. E Fo^’°(H^). ( 68 ) 

Note also that as / E L°°(M^) has compact support we then have / E Therefore we have 

existence and uniqueness of a solution U E L7q’'^(H®) to (IHHp . and estimate (IHHp follows directly from Lax- 
Milgram Lemma and from (I65p . 

Step 2. We will use [l2l Theorem 0.5]. Let rj E V(R^), r/ E [0, 1], r/ = 1 in Br and r/ = 0 in M.^\B 2 r. Then 
r]U satisfies 

-A{r]U) = fv- 2VU • Vrj - UAr] =: g in := 9.^ n B 2 R, 
rjU = 0 on dkl^R. 

From i) above g E L^(H|^). By local regularity of —A, as f ^ L°°{R^) we get U E W^’P{Br^ 2 r) for all 
p E ( 1 , 00 ), where Br^ 2 R ■= { 2 ;, ii < |x| < 2R}. Therefore, as p = 1 in it follows Vr/ = 0 in so we 
get g E W~^’P{92r) and 

\\9\\w-^’P{n^y < C{9,R,p)\\f\\L<^. 

We recall m Theorem 0.5], which states that there exists pi depending on LI^r, satisfying pi > 3 for 
every N > 3 and pi > 4 if iV = 2, such that the map U E 14))’^(f^2ij) —AU E W~^'P{92r) defines an 
isomorphism for every p E {p'i,pi). It implies that \\r]U\\iYi,pye ) < C{9, R,p)\\f\\Loa, which implies (f67p . 
□ 
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Remark 3.17 In connection with Remark 13.41 related to the interior problem, note that with our choice of 
L = L* = —A, we have here pi = =Pi- Recall that p^ > 4 if N = 2 and Pi > 3 if N > 3. 


3.4.2 Estimates of the shape derivatives of the solution 

We have the following theorem, which is similar to Theorem 13.91 


Theorem 3.18 Let be as in m- Assume it is also Lipschitz and simply connected. Let 9 £ Q, fig = 
(I + 0)(ng), / G L°°{M.^) with supp(/) C Br, Ug G the solution of 07]) (see i) in Proposition 


frm) . U 0 = Ueo{I + e). Then 

i) If f £ m £ N*, then 9 tjg £ Hq{LIq) is of class C™ in a neighborhood of 9 = 0 £ Q. 

In the following, we assume the previous hypotheses are satisfied for m = 2. 


a) Furthermore, ifLlo is Lipschitz, then for all p £ {2,Pi) where Pi is introduced in Remark\3A% we have 




[ Vf/o 

■Uo 

Jog 


/ Vt/o 

u'^ 

Jog 



^ C'l||^||^yl,2p/(p-2)(Qe^), 

< C-2||^|| 


where U( = U'g{m), K Ci 


Ci{Qo,Uo,f,R,p). 


(69) 

(70) 

(71) 


Proof. 

Step 1. The proof of i) is similar to the one of the same result in the interior, see i) in Theorem 13.91 For 
convenience we present the proof. 

The proof is based on the implicit function theorem. After changing the variable x = (I + 9){y), the weak 
form (16811 in fig is transformed in fig as follows 


f VUg ■ Mg ■ Vif - fgif = 0, 
Jni 


(72) 


where Mg = (T AT)det[Id + 'V9], T = [Id + Vd] ^ and / = / o (/ + ^(t))det[I + V0]. Consider F defined by 


F : 


0 X 

i9,U) 


= fno ■ Mg-Vp- fgp 


It is easy to check that F is well defined and of class C™ in a neighborhood of (0, Uq) G 0 x i7Q’’^(flg) (here 
we use the fact that / is with compact support), and 


duF{f), Uo){U) [ (yU ■ Vp)dx, Vt/, p £ i7o’°(f![J). 

J Qq 

Note that duF{0, Uq) defines an isomorphism from i7Q’*^(flg) to 77”^’‘^(flg), see Proposition 13.16l Then, from 
implicit function Theorem, there exists a C'” map 9 U{9) such that F{9, U{9)) = 0 for ||0||e small. Prom 

the uniqueness of solution to (|lip we obtain Ug = U{9), which proves the differentiability of 0 i—?• Ug. 
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Step 2. Again the proofs of (f69j) - (f7T|) are similar to the proofs of (except that the choice of suitable 

test functions is easier here). For the proof of (f69]l we differentiate d72|) and get 

/ {VU's-Me-Vif) = [ f'sip-VUe-M's-Vip. (73) 

With 0 = 0 and (p = as Supp(^) C B{0,R), from ([731) we get 

l|Vt>^|li2(og) < Cl / (l + |Vf/o|)(|l| + |Ve|)(|f7'| + |V17'|) 

J 

< (71(1 + ||t/o||vi/l.P(0|j))||?||iyl,2p/(p-2)(ll|^)||Vf7o||i2(Qe), (74) 

with Cl = Ci{Qo, R, f,p), which proves ([69]) . 

Similarly, if we take 0 = 0 and p = Uq in ([75|1 we get 


/ VUh-VUo < Cl (l + |V[/o|)(|e|+ |Ve|)|Vt/o| 

< Ciil + \\Uo\\ 

with Cl = Ci{Qo, f, R,p), which proves ([70]) . 

Step 3. For the proof of (f7T]l we differentiate ([7^ at 0 = 0 and take p = Uq. It gives 


Vt/^' • Vt/o 


< 


C 2 [ (1 + |Vt/o|2)(|C|2 + |vep) + |V[/o||Ve||Vf7^| 
J^R 

< (72(^||Co||wi.p(n|j)l|V.^llL2p/(p-2)(Q|^)||Vf7o||i2(fi|^) + 

(1 + llColl 

< (72(1 + ||t7o||^l,p(Qe^)||C||^l,2p/(p-2)(ll|^)) 
with C 2 = C 2 {^o, f, R,p), which completes the proof. 


□ 


3.4.3 Estimate of shape derivatives of the energy 

Theorem 3.19 Let Ll be as in m- Assume it is also Lipschitz and simply connected. Let 0 G 0, = 
(/ + 0)(n3); / £ VF™'’°°(M^), m € N*, with supp(/) d Br, Uq € the solution of (E3) and 

£{6) = \VUg\‘^. Then £ is of class near 0 = 0 G 0 and for m > 2 and for all p G (2,Pi) with Pi 

defined in Remark \3.11\ we have 

1^ (0)(^)| < (7l||^||iyl.p/(p-2)(Qe^), (75) 

i^"(o)(e,e)i < C2ii^f^ l'2p/(p-2)(fl|^)) (76) 

where Ci = Ci{Qo,Uo, f, R,p). 

Proof. The proof of the theorem is very similar to the one of Theorem 13.91 We present it briefly. Note that 

£{e) = [ VUe-Me-Vilg. (77) 

Jni 

From Theorem 13.181 it follows that £ is (7"^ in a neighborhood of 0 = 0. 
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Differentiating (f77|) at 0 = 0 in a direction ^ G 0, and then using (170]) gives 


l^'(o)(OI < 
< 


Cl 



VUo • VLf' 




+ [ |VC/oP|VC| 


with Cl = Ci{Clo,Uo, R,p, f), which proves ([7^ . 

Differentiating twice (177)1 at 0 = 0 in the direction and then using (l69]l and (f7T]l gives 


l^^"(o)(C,OI < 
< 


C2 



• VUo 


+ [ |V[>^P + |VC/oP(|C|' + |VeP) 


csiieii 


2 

tyl.2p/(p-2)(Qe)) 


( 78 ) 


(79) 


with C 2 = € 2 (^ 0 , Uo, R,p, f), which completes the proof. 


□ 


4 Application to optimal convex shapes 

In this section, we remind the strategy from |19] . and emphasize an application of the second order shape 
derivatives estimates obtained in Section E21 first in dimension two, then in higher dimensions. 

4.1 Application in the planar case 

Here is the main result of this subsection: 

Theorem 4.1 Let LIq an open convex subset ofM?, solution of the optimization problem: 

min ^j{n) = R{E{n),\n\) - p{n), n C Mfopen, convex and such that dLi C {x, a < |x| < 6}}, (80) 

where ii : —)• M zs smooth, {a,b) G (0,oo]^, a < b, and E is an energy like corresponding to the 

interior problem, or like m, corresponding to the exterior problem. 

Then every connected component o/((9Do)m := i9Don{x, a < |x| < b} is a finite union of straight segments. 

We insist here on the fact that the existence of an optimal shape for problem (I80p is true and easy to obtain 
in the case 0 < a < b < 00 , see Remark Ml The cases b = 00 and/or a = 0 require more attention: it 
may happens that minimizing sequences are not bounded (for example if ii = 0), or that they converge to 
a segment. 

Remark 4.2 This result was obtained in m for the two following particular cases: 

• E = Ef is the Dirichlet energy associated to / G R/q^(M^): 

£(0) =/^ (iiVCInP -/%) , where [/„ eolves { 

• £' = Ai is the first Dirichlet eigenvalue of D: 

Ai(D)=min|£lJ^, [7 G Ro'(^^) \ w} • 
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Remark 4.3 In Theorem 14.11 we chose a simple function of P{Q). Actually, the same result would extend 
to functionals of the form R(^E{Q), |n|, Ai(n), P{^)) with p —)• R{-, •, ■,p) being decreasing and concave. 


Remark 4.4 Similarly to |191 Theorem 4, Example 8], Theorem 14. II is also valid with a volume constraint. 
More precisely, any solution of 

min{J(n) = R{E{Q,)) — P{^), C open, convex and such that |n| = Vq} (81) 

(where R : M —)• M is smooth and E is as in Theorem 14.ip is a polygon. See also Remark 14.61 for other 
geometrical constraints. 

In order to analyze problem (|80l) . we use, as in |18l 119] . the following classical parametrization of 2- 
dimensional convex with polar coordinates {r,9) E [0,oo) x T, where T = M/27rZ: 


Hn := \ {f,9) E [0,oo) X R ; r < 


1 


u{e) / ’ 

where u is a positive and 27r-periodic function. A simple computation shows that the curvature of is 


(82) 


i^anu = 


u" + u 




(83) 


This implies that Qu is convex if and only if u" + u > 0, which has to be understood in the sense of H ^(T) 
if u is not More precisely, if u E R^(T) then u" + u > 0 if and only if 

V u E R^(T) with u > 0, / {uv — u'v') dd > 0. 

JT 

Throughout this section, any function defined on T is considered as the restriction to T of a 27r-periodic 
function on R, with the same regularity. 

With this parametrization, considering j{u) = J{Qu), Problem (|80p is equivalent to 

min {j{u), u" + u > 0, u E Uad} > where Uad = {w E IT^’°°(T), l/u E [a, b]}. (84) 

Then we have the following result proven in m Theorem 3]. 


Theorem 4.5 Letuo > 0 be a solution for (l84ll andTin := < 0 E T, a < —Assume j : IT^’°°(T) —>• 

I M&) } 

R is and that there exist s E [0,1), a > 0, /3,7 E R such that, for any v E IT^’°°(T), we have 

f'{uo){v,v) < +7bl//i(T)ll^^llHRT) +/3 ||^^|Ihs(t). (85) 

If I is a connected component o/Tj„, then 

u'o + uq is a finite sum of Dirac masses in I. 

This result, combined with the estimates from Section 12.21 will lead to a proof of Theorem 14.11 Indeed, 
formula (|83l) explains that Du is polygonal if and only if u" + n is a sum of Dirac masses. 


22 





Remark 4.6 As in m, we may focus on a geometrical constraint different from dil, C {x,a < \x\ < b}. In 
that case, the previous results remain valid when replacing the definitions of (9no)m with: 

Tin = {6* E T / Be, 5 > 0,Vn E — £,0 + e) with ||u||yi/i,oo < <5, uq + u E Uad}- 

E Srio / 30 E Tin,x = —^(cos0,sin0)| . 

where Uad replaces {n/9n C {x,a < \x\ < 5} and includes the geometrical constraint, except the convexity 
one. Note that (df}o)in describes the part of the boundary which is inside the constraints, apart from 
the convexity one. However, if one deals also with a finite-dimensional equality constraint (for example 
1^1 = To); then one should work with a Lagrange multiplier, see Remark 14.41 



Remark 4.7 Due to the convexity constraint, the existence of a solution to ()80p can be proved easily when 
0 < a < b < oo. We briefly provide a proof here: let (Hn)nGN £ parametrization E Uad, be 

a minimizing sequence of J(H). From u'^ + Un > 0 and l/un E [a,b], it follows that {un)n£N, is strongly 
relatively compact in H^(T) and therefore we may assume linin^ooUn = no in H^{T) (see also Remark 
I4.12|] . with Uq B-uq >0 and I/uq E [a, b]. In such conditions the domains satisfy a uniform exterior cone 
condition and will converge to Hq for the Hausdorff convergence, with parametrization uq. Then will 
converge to I/qq in iL^(H(0, 6)) (see [!!( Prop 2.4.4 and Theorem 3.2.13]), and thanks to the continuity of 
the perimeter for the H^(T)-norm we get that P(H„) converges to P(Ho)) which implies that Hq solves (|80ll . 

Proof of Theorem 14.11 Let no such that Hq = Then uq is a solution of (I84p . where j(n) = 

R{e{u),m{u)) — p{u) and e(n) = E{Q,u),^{u) = \idu\,p{u) = P(H„). 

The second order derivatives related to the geometrical terms m(n) and p{u) can be easily computed. 
Indeed, first note that for every u E IF^’°°(T), rt > 0 we have 

Vu'^ + n'2 

-5-40. 

n^ 

It follows that m and p are twice differentiable in VF^’°°(T) and there exist some real numbers j3\, 1^2,1 
(depending on uq) and a > 0 such that, 


= L a?'®’' 


p{u) 


\m''{uo){v,v)\ < /?i||n|||2(T) 

p''{uo){v,v) > - i\v\h^(t)\\v\\l2(j) - P2\\v\\l2(^j^ 


It remains the difficult term e”{ uq){v,v), which relies on Section 12.21 To handle this term, we consider 
V E iy^’°°(T) and introduce the vector held: 


^EC2((-to,to),lTb“(M^M2^), 

^(t) = ^—-j—^- ^e*^77 (r, 0) on (in polar coordinates). 


( 86 ) 


where p E C“(M^) = D(M^),0 < p < 1, p = 1 (resp. 77 = 0) in a neighborhood of SHq (resp. the origin). 
Let us point out that (/ + ^(t))(H„) = Q-u+tv, e(n + tv) = E{Q.u+tv) = £{i{t)) and also that 


2 

Mv E ITb°°(T) : e'(0)(w) = ,e"(0)(w,?;) = on 5H. 

Then, we will differentiate twice around 0 the map 1 ip{t) := E{i}u-\-tv) = 

^'( 0 ) =e'{uo)iv) =£'mem, 

p^'iO) =e"{uo){v,v) =£:"(0)(e'(0),^'(0))+f'(0)(e"(0)). 


( 87 ) 
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and therefore Corollary 12.61 (see also ([2^ . (l23]l l implies there exists e E (0,1/2) such that 

|e"(no)(^,z;)| < +/34(£)||e'(0)||^l-e(ao)■ 

We then use the fact that the and norms are equivalent (since the transformation 

V’ = ip{r, 9) := is bi-Lipschitz near T and V’(l') = dO), and the fact that is a Banach algebra 

to obtain 

\e"{u){v,v)\ < /35(e)IIv||^l-e(T), 

Easy computations about the term /^(^(n), |fl|) imply that j" satisfies ([851) . and therefore Theorem 14.51 
applies. It follows that Uq + uq is a finite sum of Dirac masses in any connected component of Tj„. Geomet¬ 
rically speaking, considering the formula (|83p . this correspond to the fact that any connected component of 
((9Do)m is polygonal and concludes the proof. □ 

4.2 Application in the multi-dimensional case 

In the multi-dimensional case, convexity constraint in shape optimization is much less understood, see 
mm and the work of T. Lachand-Robert. Nevertheless, we can use a parametrization similar to the one 
used in Section liTl and we show in this section that our estimates of shape derivatives allows to obtain 
results in any dimension. 

For A > 2, if u : —)• (0, 00) is given, = {x E |x| = 1}, we can consider 

:= |(d6 ') € [0,oo) X 

The function u is the so-called gauge function of The set is convex if and only if the 1-homogeneous 
extension of u, denoted by the same letter and given by u{x) = |x|u(x/|x|) is convex in (in this section, 
we will refer to this property by saying that u : —)• R is convex), see m Section 1.7] for example. 

In this way, we describe every bounded convex open set containing the origin. Throughout this section, 
the regularity of any function defined on is seen as the regularity on R'^ \ {0} of its 1-homogeneous 

extension, and it is classical that it is equivalent to the regularity of the set itself. 

With this parametrization, considering j{u) = J(D„), problem ([T]) is equivalent to 

min {j{u), u : —)• (0, 00 ) convex , u E Uad] , (89) 

where Uad C IT^’°°(S^“^) is defined as Uad = {u E IT^’°°(S^“^), E Oad} and Oad is the class of 

admissible open bounded sets considered in ([T]) (not taking into account the convexity constraint). 

Then in the same spirit as Theorem 14.51 we can prove: 

Theorem 4.8 Let uq > 0 be a solution for (f89|) . where Uad is a convex subset of IE^’°°(S^“^). Assume 
j ; VE^’°°(S^“^) — R is and that there exists s E [0,1), a > 0, /3 ,7 E R such that for any v E IT^’°°(S^“^) 
we have 

j”{uo){v,v) < -a|u|^i(gjv-i) +7klj^i(§iv-i)||u||j:/s(§iv-i) -f/3||u||^,(giv-i). (90) 

Then the set 

Tuq = {u E IT^’°°(S'^“^)/3e > 0, V|t| < e, uq + tv ^ Uad o.nd is convex}, (91) 

is a linear vector space of finite dimension. 
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Proof. It is easy to check that Tug is a linear vector space. Since uq is optimal, we have j"{uo){v, v) > 0 for 
every v G To conclude, we prove that the unit ball of for the is relatively compact. 

Indeed, thanks to ([90]), we get 

Vn G Tuo such that ||u||j|^s(§iv-i) < 1, Q;|u|^i(gjv-i) < + /?, 

which easily leads to 


Vu G Tuq such that ^ 1) 




ItI + a/7^ + 4«|/3| 
2a 


From the compact embedding of in we conclude that the unit ball of C H^{S^ 

has a compact closure and therefore from Riesz Theorem, T^q has finite dimension. □ 

As a corollary, the results we obtained in this paper, namely iii), Theorem 12..SI lead to the following 
generalization and improvement of [S] Theorem 4.5]: 


Corollary 4.9 Let Qq be an open eonvex subset solution of the optimization problem: 

min = R{E{n),\n\) - p{n), open, convex and sueh that dfl C {x,a < jxj < 6}|, (92) 

where A : —)• M is smooth, {a,b) G (0,oo]^, a < b, and E is an energy like |3), corresponding to the 

interior problem. 

Then denoting by uq the gauge function of Qq, defined in (19X1) has finite dimension. In partieular, 
if io is a relatively open subset of {dQo)in := dLlo (7 {x,a < |x| < b}, then the Gauss curvature of Qq 
vanishes on u. 

Compared to [3l Theorem 4.5], we have enlarged the class of functionals under consideration, and we require 
less regularity on the optimal set. 


Remark 4.10 Contrary to Theorem 14.11 we cannot deal with the exterior problem when N > 2. This is 
due to the fact that for a Lipschitz domain LI in dimension N > 3, our estimates do not imply that there 
exists s < 1 such that (whereas it is the case for N = 2, see (l23]) l. 

Remark 4.11 As in Eemark 14.61 we may focus on a geometrical constraint different from dLl C {x,a < 
|x| < b}. In that case, given Oad a class of admissible open bounded set, the previous result remain valid 
when replacing the definition of (cAlo)m as follows 

6*1 , with 

= {6* G / Be,6 > 0,Vu G WQ’°°{B^N-i{0,E))with IJuUvi/i.oo < 6,uo + v & Uad} and 

Had — € Glad}i 

if Had is convex, as assumed in Theorem 14.81 


{dLlo)in = < X G dLlo / 36* G S, 


N-l 
in ' 


X = 


uo{9) 


Remark 4.12 As in Remark 14.71 the existence of a solution to ()92D can be proved easily when 0 < a < 6 < 
oo. The proof is similar to the one in Remark 14.71 the only adaptation to the multi-dimensional case is the 
fact that a sequence of 1-homogeneous functions : Ri —)• M (where Bi is the unit ball of M^) such that 
Vn G N,V6* G < Un{0) < (5, with 0 < a < /3 < oo, is strongly relatively compact in H^{Bi) (actually 

in VF^’^(Ri) for any p < oo). To that end, we first notice that IJVunllcxD < /3 (see for example 0), which 
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implies that up to a subsequence Un converges to u for the L°°(i?i)-norm. Clearly u is then 1-homogeneous, 
convex, and satisfies MO G < u{6) < /3. To obtain the convergence for the //^-norm, we follow [21 

Lemma 2.1.2]. As in Remark l4.71 the perimeter, the volume, and the energy E are continuous, and therefore 
the set whose gauge function is the limit of Un is a minimizer. 

Proof of Corollary 14.91 Let uq such that flo = ^uo- Then uq is a solution of (j89l) . where j{u) = 
R{e{u),m{u)) — p{u) and e{u) = E{Qu),m{u) = \Qu\,p{u) = P{^u)- 

The second order derivatives related to the geometrical terms m{u) and p{u) can be easily computed. 
Indeed, first note that for every u G iy^’°°(S^“^), u > 0 we have (see for instance m) 


m{u) = / 


1 


-1 Nu^ 


dO, 


p{u) = / 

Js^-i 


+ \VrU\^ 


dO. 


where V,- denotes the tangential gradient on It follows that m and p are twice differentiable and 

there exists /3i (depending on uq) such that. 


\m”{uo){v,v)\ < /3i||u|||2(t) 

Easy direct calculations give 

p{u){v) = f duG{u,\Vru\)v + dpG{u,\Vru\) 


VtU ■ VtV 


p''{u){v,v) = 


[ duuGi 


VrU\ ’ 
VrU • Vt-V 


u,\Vru\)v^+ 2dupG{u,\VrU\)- 

VrU 


-V + 


O IV7 -VrU)^ ^ ^ ^ - (V^U • V.rU)2 


Taking into account dpG{u,p) = 
terms gives 


P 


iN 


Vrltj 
, dppG{u,p) = 


{v? + p2)3/2 


|V^«|3 

and rearranging the last two 


p''{u){v,v) = f duuG{u,\Vru\)v‘^ + 2dupG{u,\VrU\) 

Js^-i 


VrU • Vt-U 


Vt-U 


-V + 


+ jVT-ttj 


iV.uL - 


|V.rril 


+ I Vt-uP V I Vt-U 


Vt-u 


VtV 


As the first and second terms are controlled by /32||u||^2 and 7 ||u||j;^ 2 |u|j:^i respectively, with /32,7 depending 
on uq, and 


IV7 |2 

V.r'f — 


VtU ) > (1 — (1 — a))|Vrr’|^ = a||VT-f 1^, 


Vr^ 

u'^ + |VT-rtp \ IVt-w] 

with a = a{uo) > 0, it follows 

p"{uo){v,v) > a|v|^i(§iv-i) -7l^lHi(s^-i)ll^lli2(s^-i) -/32||?^|li2(giv-i). 

As in the proof of Theorem 14.11 the term e"{uo){v,v) relies on Section [2.21 this time in any dimension, but 
restricted to the case of convex domains. With the same proof and using iii). Theorem 12.31 we obtain 

|e"(.o)(.,.)| < /33WII^"(0)llH..-L7ao) +/54(r)||?'(0)||;,. . 
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valid for any r G (1, oo). We would like to write the above inequality in terms of the ff*(§'^“^)-norm of v, 
for some s G (0,1). For this, first we note that any bi-Lipschitz transformation from to Sflo (in our 

case it is 0 {l/u{6),9)) defines a diffeomorphism from W^’P{dQo) to for every s G [0,1] and 

p G [l,oo]. Then the above inequality is equivalent to 


e''{uo){v,v)\ < I35{r,u) 


Un 




+ Pair, u) 

,r-(gJV-l) 


V 



2 


(93) 


Note that here we used the formulas ()87p with e*® replaced by 6 (of norm 1). 

For the first term in ()93p we wonder for which s G (0,1) the following inclusion is continuous: 


H 


sf-mN-l 




(94) 


We use [26l Corollary, p. 189], which refers to the more general embedding C ^ in . 

Note that by using a partition of unity, some Lipschitz transformations and the extension theorem, we may 
use these embeddings as if was and therefore we use it with N — 1 instead of N. In our specific 

case, this leads to the condition 



N 

Y 


< s < 1. 


(95) 


Clearly, for every s G (0,1) one can find r so that (|95p holds, which implies that (j94p is valid. 
For the second term in (|93p . we look for the inclusion 


(96) 


which is valid (again using [26]) if 

l-:^ + ^^(l--)<^<l- (97) 

2r 2 r 

Inequality (fWp has a solution r (sufficiently close to 1) for any s G (1/2,1), which implies that (fMp holds 
for every s G (1/2,1). 

Combining (f^ and (IMp with ([Mil gives 

\e"{u){v,v)\ < /37(s,Mo)lbllH«(siv-i), 
for any s G (1/2,1). 

Easy computations about the term B(B(n),|n|) implies then that j"{uo) satisfies (IMIl . and therefore 
Theorem 14.81 applies, and Tuq has finite dimension. 

The fact that the Gauss curvature must vanish where it is defined is an easy consequence. Indeed, if u 
is as in the statement of Corollary 14.91 and if the Gauss curvature is positive at some point of lo, then it is 
positive in a neighborhood a) of this point. As a consequence, any smooth function with compact support 
in Lu is in Tug- This contredicts the fact that has finite dimension. □ 
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